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Abstract 


2 Episodic events of both Saharan dust outbreaks and African Easterly Waves (AEWs) are 

3 observed to move westward over the eastern tropical Atlantic Ocean. The relationship between 

4 the warm, dry, and dusty Saharan Air Layer (SAL) on the nearby storms has been the subject of 

5 considerable debate. In this study, the Weather Research and Forecasting (WRF) model is used 

6 to investigate the radiative effect of dust on the development of AEWs during August and 

7 September, the months of maximum tropical cyclone activity, in years 2003-2007. The 

8 simulations show that dust radiative forcing enhances the convective instability of the 

9 environment. As a result, most AEWs intensify in the presence of a dust layer. The Lorenz 

10 energy cycle analysis reveals that the dust radiative forcing enhances the condensational heating, 

1 1 which elevates the zonal and eddy available potential energy. In turn, available potential energy 

12 is effectively converted to eddy kinetic energy, in which local convective overturning plays the 

13 primary role. The magnitude of the intensification effect depends on the initial environmental 

14 conditions, including moisture, baroclinity, and the depth of the boundary layer. We conclude 

15 that dust radiative forcing, albeit small, serves as a catalyst to promote local convection that 

16 facilitates AEW development. 
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1 1. Introduction 

2 Over the eastern tropical North Atlantic Ocean, episodic outbreaks of mineral dust off the 

3 African continent are observed to move westward in association with the propagation of African 

4 easterly waves (AEWs). The dust plume and the AEWs straddle the African easterly jet (AEJ) 

5 about 15°N (e.g., Karyampudi and Carlson, 1988; Karyampudi et al., 1999). The dust plume, 

6 which is embedded in the Saharan air layer (SAL) — the westward moving warm parcel of air 

7 that develops over the African continent — generally extends between 10°N and 25°N, while 

8 AEWs propagate south of the AEJ, between 5°N and 15 °N. The AEWs usually have a period of 

9 3-8 days, wavelength of 2000-4000 km, and phase speeds 6-8 m s' 1 (e.g., Reed et al., 1977; 

10 Karyampudi and Carlson, 1988; Karyampudi et al., 1999; Zipser et al., 2009). Satellite images 

1 1 show convective clouds around the trough of the waves indicating the convectively active part of 

12 the AEW in association with positive relative vorticity at 700 hPa (e.g., Berry et al., 2007). The 

13 mechanism for the development of AEW was originally explained by the barotropic-baroclinic 

14 instability theory (Charney and Stem, 1962; Burpee, 1971). In recent years, Hsieh and Cook 

15 (2007; 2008) found that AEWs could be developed through the baroclinic overturning over 

16 ocean in the absence of shear instability associated with the AEJ. When environmental 

17 conditions are favorable for development, an AEW can intensify and eventually grow into a 

18 tropical cyclone (TC). It has been documented that most major hurricanes (> category 3) are 

19 originated from AEWs (Landsea, 1993). 

20 The fact that the dust plume is often in the vicinity and at similar altitude of the vorticity 

21 center of the AEWs has generated a lot of discussions on the linkages between the SAL and the 

22 AEWs. For example, using geostationary satellite data, Dunion and Velden (2004) found that 

23 Hurricane Joyce weakened significantly right after the SAL reached Joyce while Hurricane Issac 
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1 curved northwest and avoided the SAL, and intensified into a category-4 hurricane. Evan et al. 

2 (2006) established an inverse correlation for the years 1982-2005 between atmospheric dust 

3 cover and TC activity as measured by satellite. Lau and Kim (2007) also reported an inverse 

4 correlation between the aerosol index from Ozone Monitoring Instrument and the Sea Surface 

5 Temperature (SST), suggesting that the cooler sea surface due to reduced insolation by dust 

6 could suppress TC activity. Reale et al. (2009) found that the atmospheric thermal structure 

7 associated with the SAL is unfavorable for the AEW development. 

8 Karyampudi and Pierce (2002), on the other hand, found that the SAL had a positive 

9 impact on two of three hurricanes examined, and a negative impact on one case which was in a 

10 particularly dry year. They proposed that the SAL could increase the meridional temperature 

1 1 gradient, and, thus, enhance the baroclinic instability on the leading and southern edge of the 

12 SAL. The baroclinic instability can be released to facilitate further development of the 

13 disturbance. Koren et al. (2005) analyzed Moderate Resolution Imaging Spectroradio meter 

14 (MODIS) data and reached a similar finding, that the deep convection and aerosols are positively 

15 correlated. Moreover, Jones et al. (2004) found that the lower troposphere associated with the 

16 AEWs is warmed, which favors AEW intensification. They hypothesized that the dust radiative 

17 effect is responsible for such warming, corresponding to the findings in Alpert et al. (1998). 

18 Nevertheless, Braun (2010) analyzed multiple satellite datasets and found that the dust plume 

19 rarely reached the location of the disturbance to affect the AEW development. He found that the 

20 dry air responsible for suppressing TC development was a result of large-scale subsidence 

21 associated with semi-permanent high-pressure systems, and the SAL did not play a key role in 

22 this context. 
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1 Although previous observational studies have not agreed on the influence of the SAL on 

2 AEWs, numerical models have been implemented to investigate the effect of the SAL. One 

3 advantage of the modeling approach is that one can identify the effect of one factor/process of a 

4 complex system while keeping all other factors unperturbed. Many modeling studies have 

5 focused on the effect of aerosols in the SAL. The microphysical effect of aerosols on TCs has 

6 been suggested to have a negative impact on TC, but the mechanism involved is not yet 

7 conclusive. Studies show that increasing aerosols in the environment can weaken a TC by 

8 increasing the vertical wind shear (Fan et al., 2009), blocking the energy flow into the core of a 

9 storm by enhancing convection in the outer rainbands (Zhang et al., 2009; Carrio and Cotton, 

10 201 1), and diminishing convection in the core (Khain et al., 2008). 

1 1 The direct radiative effect of aerosols has been found to influence the AEJ and 

12 convection in the region of interest (e.g., Tompkins et al., 2005), but its role in AEW 

13 development is not entirely clear. Lau et al. (2009) used a global GCM and found that the dust 

14 radiative effect could fashion an “elevated heat pump" mechanism (Lau and Kim, 2006) to cool 

15 the surface, warm the mid-troposphere, and create a circulation pattern that enhances subsidence 

16 over TC’s main development region in the tropical North Atlantic. All of these effects are 

17 unfavorable for the development of TCs. Sun et al. (2009) found that the radiative effect of the 

18 SAL contributes to an atmospheric thermal structure that is unfavorable for AEW development. 

19 These studies are consistent with Dunion and Velden (2004), Evan el al. (2006), Lau and Kim 

20 (2007), and Reale et al. (2009). In Stephens et al. (2004), the authors found that deep 

21 convections are prohibited in the first 10 days of simulations due to dust radiative forcing. 

22 However, increased occurrence of deep convection in the dust -covered region compared with the 

23 dust-free region is observed after 10 days of simulation, which supports the observational 
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1 analysis in Koren et al. (2005), and is consistent with Karyampudi and Pierce (2002) and Jones et 

2 al. (2004). 

3 In this study, our focus is on the direct radiative effect of aerosols on the development of 

4 AEWs, regardless of whether or not the disturbance eventually grows into a TC. Using the 

5 Weather Research and Forecasting (WRF) model, numerical experiments were conducted for 60 

6 observed AEWs, all occurrences that lasted at least 4 days within the domain (5°N-25°N; 10°W- 

7 50°W) in August and September, in years 2003-2007. For each of the AEW cases, we conducted 

8 three simulations, one without the radiative effect of the dust, one with a dust layer located 

9 between 700 hPa and 850 hPa, and the other one with a dust layer located between 600 hPa and 

10 850 hPa. A Lorenz energy cycle analysis (Lorenz 1955; 1967) was then applied to understand 

1 1 how aerosol direct radiative forcing modifies the energetics of the atmosphere and, in turn, 

12 influences the development of the AEW. In this paper, the synoptic conditions are described in 

13 Section 2. A brief description of the model is given in Section 3. Results from the numerical 

14 simulations of 60 AEWs are presented and discussed in Section 4, and conclusions are drawn in 

15 Section 5. 

16 

17 2. Synoptic conditions 

18 To establish the climatology of the AEW but to avoid effects from large-scale climatic 

19 oscillation such as ENSO (El Nino-Sothem Oscillation), we limit our study to only 5 consecutive 

20 hurricane seasons, from 2003 to 2007. MODIS Aerosol Optical Thickness (AOT) was obtained 

21 from the Level 1 and Atmosphere Archive and Distribution System 

22 (http://ladsweb.nascom.nasa.gov/) at NASA/GSFC. It is Level 3 processed data, consisting of 

23 the daily average on a 1° x 1° latitude-longitude grid of the 10 km x 10 km retrievals from the 
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1 MODIS sensor on the Terra satellite (Remer et al., 2005). Relative vorticty, winds, and 

2 temperature were obtained from the ECMWF operational analysis for the same period. For 

3 comparison with MODIS observations, we use only the 12Z ECMWF data, in order to maximize 

4 alignment with the Terra satellite, which passes over the eastern tropical North Atlantic Ocean in 

5 the late morning. 

6 The data shows that in the boreal summer, the strongest easterly wind, which defines the 

7 AEJ, is observed at about 15°N (Figure 1, left panel). In the meantime, positive relative 

8 vorticity, which indicates the vorticity field of the background flow as well as the trough of the 

9 AEWs, are in the area south of the AEJ and south of the dust plume, which is north of the AEJ. 

10 Relative vorticity is strongest near the coast of Africa, and its magnitude gradually decreases as 

1 1 the AEWs propagate westward. The dust layer extends appreciably from 10°N to 25°N and 

12 spreads up to around 40 degrees west of the African coast, with a maximum close to the coast. 

13 As the dust plume travels westward, AOT decreases due to settling on the ocean surface. 

14 Regarding inter-annual variability, the dust plume is the most opaque and extends the farthest in 

15 2003, and the least opaque in 2007. 

16 The AEWs poleward of 15°N are dominated by the relative vorticity field at 850 hPa and 

17 AEWs equatorward of 15°N are dominated by 600 hPa relative vorticity (Thorncroft and 

18 Hodges, 2001). Following the work by Berry et al. (2007), we chose 700 hPa to analyze the 

19 AEWs in this study. The temperature and winds at 700 hPa (Figure 1, middle panel) shows that 

20 the warm air to the north and the cool air to the south, straddling the AEJ at about 15°N, result in 

21 horizontal wind shear about the AEJ. The anti-cyclonic circulation associated with the warm air 

22 parcel indicates a convectively stable condition with subsidence in the area. Due to the thermal 

23 wind relation, the strongest easterlies are found at about 700 hPa. The warm air extends to about 
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1 40°W in most years, except in 2004, when the warm air was confined mostly east of 20°W. The 

2 meridional temperature gradient is about 4 K over 20 degrees of latitude over land and only 

3 about 1 K over ocean at 30 °W. At about 40°W, where the environmental shear instability 

4 vanishes, baroclinic overturning could reverse the sign of the gradient of meridional potential 

5 vorticity and thus provide instability (Hsieh and Cook, 2007; 2008). The meridional temperature 

6 gradient at 1000 hPa (Figure 1, right panel) varies from 16 K per 20 degrees over land (at 5 °W) 

7 to -2 K over ocean, where the air temperature near the surface is strongly affected by the SST. 

8 The data shows that the warm air associated with the SAL is elevated in altitude and does not 

9 extend to the surface. 

10 In Figure 2, the time series of 700 hPa relative vorticity and AOT, averaged between 

1 1 15°W and 50°W, shows that the strong relative vorticity to the south is often accompanied with 

12 high AOT to the north. The correlation coefficient between the time series of AOT averaged 

13 over the north part of the domain and the time series of relative vorticity averaged over the south 

14 part of the domain is 0.34. TC-genesis often occurs at the southern edge of the dust plume and 

15 the northern edge of the vorticity field of the background flow. Furthermore, The amplitude and 

16 latitudinal positions of maximum AOT and relative vorticity vary with season in tandem. In 

17 June and October, the dust plume and the AEW are in further south than in July, August, and 

18 September. Both dust and AEW activity are generally calm in October, when TC-genesis is less 

19 frequent. The dust outbreaks are stronger and more frequent in the early months than in the latter 

20 months of a hurricane season, except for 2005, in which the dust outbreak events took place 

21 mostly from July to mid-September. These features are also discussed in previous studies (e.g., 

22 Calrson and Prospero, 1972; Dunion, 201 1). On the other hand, the AEW tends to be weaker in 


23 June, July, and October, and stronger in August and September. Although one would expect a 



1 more intense disturbance to have a better chance of becoming a TC, not every intense AEW 

2 disturbance reaches that stage. More TCs are developed in 2003, 2004, and 2005 than in 2006 

3 and 2007, despite the fact that the number of AEWs and dust outbreaks did not change much 

4 (Table 1). 

5 Since most TCs occur in late summer, we limited the period for the numerical 

6 simulations to August and September, and identified a total of 60 AEW disturbances over the 5- 

7 year period. There was only one criterion for selection, viz., the AEW last at least 4 days. The 

8 700 hPa relative vorticity of the AEW must be larger than 1 x 10~ 5 s' 1 for at least four days within 

9 this domain in order to be considered in this study. The dates that the AEWs leave the coast of 

10 Africa are shown in Table 2. As one can see in Figure 3, on average, the vorticity center of the 

1 1 AEWs moves about 6 degrees westward per day. The AEWs weaken as they propagate further 

12 west. The dust plumes are generally found in the northern periphery of the vortices. Although 

13 the major part of the dust layer lies north of the vortex, it can extend as far south as the vortex 

14 center and even beyond. Figure 3 shows that part of the dust plume gradually extends from north 

15 to west of the vorticity center, counter-clockwisely in accordance with the drift of the vortices. 

16 The southern edge of the dust plumes can extend to about 8 °N, about 5 degrees further south to 

17 the composite vorticity center. The synoptic conditions of the AEWs and the associated dust 

18 layer in the 5 late summers are similar to previous observational studies and will be used as we 

19 construct our numerical experiments. 

20 It is worth noting that MODIS AOT is not derived in cloudy areas and convective clouds 

21 associated with AEWs (Berry et al., 2007) can obscure the AOT retrieval. Hence, the dust plume 

22 in the real-atmosphere can be even closer to the vorticity center (i.e., further south) because of 

23 the bias inherent in the AOT retrieval. 
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3. Model description 

The WRF Version 3.1.1 is used to investigate the direct radiative effect of dust on the 
AEW. We implemented in WRF the NASA/GSFC (National Aeronautics and Space 
Administration/Goddard Space Flight Center) radiative transfer model (Chou and Suarez, 1999; 
2001) with recent improved features (Matsui et al., 2007; as described in Shi et al., 2010). In this 
configuration, aerosol radiative properties are part of the input to the radiative transfer model, 
along with temperature and the radiative properties of the gaseous atmosphere, so that the 
aerosol radiative forcing is included in the radiative heating rate calculation. The non-radiative 
effects of aerosols (e.g., hydrological or mechanical) are not included in the model for this study 
and are beyond the scope of this paper. 

The physical parameterizations selected for the model runs are the Goddard microphysics 
scheme (Tao and Simpson, 1993; Tao et al., 2003), Kain-Fritsch cumulus scheme (Kain, 2004), 
Yonsei University planetary boundary layer (PBL) scheme (Hong et al., 2006), and Noah land- 
surface scheme (Chen and Dudhia, 2001). For initial and boundary conditions, we use the 
ECMWF (European Centre for Medium-Range Weather Forecasts) advanced operational 
analysis, obtained from the ECMWF Data Server 

(http://www.ecmwf.int/products/data/archive/index.html). The gridded analysis dataset presents 
the AEW reasonably and has been used to study the AEWs (e.g., Reed et al., 1988; Thorncroft 
and Hodges, 2001). 

It is worth noting that the warming or cooling of air by absorption and emission of 
radiation energy is instantaneous, while warming or cooling of the sea surface occurs on a much 
longer time scale due to the latency in ocean response (e.g., Arking and Ziskin, 1994). Hence, in 
our model, the sea surface temperature (SST) is not affected by direct radiative effect of aerosols. 


10 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


However, in the time scale shorter than a few days, SST is a function of wind speed. While 
strong surface winds are expected in the presence of an AEW or TC, storm- forced SST cooling 
(Price, 1981; Sanford et al., 1987) can weaken the parent (Black and Holland, 1995) and the 
future storms (Brand, 1971). In our simulations, we use ECMWF operational analysis, which 
has assimilated observational data and contains both SST and winds, as our initial and boundary 
conditions. The model calculates the meteorological fields within the domain. However, the 
SST, as the lower boundary condition, is not affected by the simulated meteorological fields. 
Therefore, the instantaneous response of SST to the high surface winds associated with the storm 
cannot be captured in our model and is beyond the scope of this study. 

The model grid is 300 (east-west) by 180 (north-south) by 24 (vertical), covering the 
eastern tropical North Atlantic Ocean where the AEW and the dust-laden SAL that develop over 
the African continent move generally westward in this region. The projection method used to 
create the domain is Lambert conformal conic projection with the two intersecting true latitudes 
being 10°N and 20°N, and the center point 15°N, 30°W. The corner points of the domain are: 
southwest 2.15°N, 49.72°W; southeast 2.15°N, 10.28°W; northwest 25.98°N, 52.07°W; 
northeast 25.98°N, 7.93°W. The horizontal grid spacing is 15 km. 

In this study, our goal is to understand the radiative effect of mineral dust on the AEWs, 
not the model’s ability simulating a dust plume. The aerosol radiative properties depend on 
aerosols’ physical and chemical properties such as size distribution, mass, composition, etc., as 
well as the conversion from aerosol’s physical and chemical properties to the radiative properties 
in the model such as refractive indices assumed and the averaging rule adopted (Fast et al., 2005; 
Barnard et al., 2010). Accurately simulating aerosol’s mass, number, and size distribution that 
involves complex processes such as the mobilization, transport, and deposition of mineral dust 
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1 remains a difficult task. As a result, the large uncertainty associated with the dust direct radiative 

2 forcing can lead to a different impact upon the climate system (Zhao et al., 2010). These 

3 uncertainties require further investigations and are beyond the scope of this study. While 

4 acknowledging the fact that there exist large uncertainties associated with the radiative property 

5 as well as the distribution of Saharan dust, we focus on the atmospheric response (i.e., the 

6 intensification of AEWs) to a prescribed dust layer, which has the properties and distribution 

7 within the range of reported dust events. 

8 For each of the 60 AEW cases (Table 2), we conduct three 6-day simulations, the control 

9 experiment of “no aerosol” (NA), a “with aerosol” (WA) experiment with a prescribed dust 

10 plume located at lower altitudes (WAl), and a WA experiment with a prescribed dust plume 

1 1 extends to higher altitudes (WAh). The dust layer is placed to the north of 15 °N, in accordance 

12 with the observation (Figure 3). One limitation of this prescribed dust layer is that the transport 

13 and removal processes of aerosols are ignored. The prescribed dust layer exerts the same 

14 radiative forcing constantly throughout the entire simulation period, resulting in a continuous 

15 forcing that gradually builds up its effect on the AEW and the environment. In the real 

16 atmosphere, mineral dust can be gradually removed from the atmosphere through dry deposition 

17 and wet scavenging as the dust plume moves further west into the Atlantic Ocean, as shown in 

18 Figure 3. As a result, its radiative forcing is expected to be the strongest when the AEW and the 

19 dust plume are close to the coast of the African continent, and fade away over time. For vertical 

20 distribution, dust aerosols are evenly distributed between 700 hPa and 850 hPa for WAl, 

21 corresponding to the high scattering layer between 1.2 km and 3 km during a vertical descent 

22 through the SAL reported in Zipser et al. (2009). For WAh experiments, the dust layer is placed 

23 between 600 hPa and 850 hPa, similar to the vertical distribution reported in Myhre et al. (2003). 
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1 The radiative properties of the dust layer are set in accordance with values in previous 

2 studies. For optical thickness, we set the shortwave AOT (r sw ) to 1.0, which is about 30% less 

3 than what Haywood et al. (2003) reported in an event with the heaviest dust loading (x sw = 1.5). 

4 The longwave AOT (r lw ) is set to 0.5, half of the x sw , as suggested in Haywood et al. (2005) 

5 while using the refractive index from Volz (1973). Note that because we prescribed the radiative 

6 properties of aerosols (optical thickness, single scattering albedo, and asymmetry factor), the 

7 refractive index is not directly used in our calculations. However, it is important to acknowledge 

8 that using different refractive indices gives very different longwave-to-shortwave (IR-to-VIS) 

9 AOT ratio. Refractive index depends on the composition of dust particles. The refractive index 

10 used in climate models is often based on World Climate Program (WCP, 1986) or Volz (1972). 

1 1 Fouquart et al. (1987) used measurements of Saharan dust over Niger for estimating refractive 

12 index. In Haywood et al. (2005), the authors reported that the IR-to-VIS AOT ratio is 0.33 while 

13 using the refractive index in Fouquart et al. (1987), and 0.52 while using the refractive index in 

14 Volz (1972). Highwood et al. (2003) stated that the as AOT at 0.55 micron is 0.67, AOT at 10 

15 micron varies from 0.4 to 0.23 when different refractive index is used, rendering the IR-to-VIS 

16 AOT ratio ranges from about 0.60 to 0.34. In Hansell et al. (2010), the authors reported IR-to- 

17 VIS AOT ratio as high as 0.75 occasionally, but adopted the refractive index in Volz (1972) for 

18 their radiative transfer calculation because their focused area is around Cape Verde of which the 

19 mineralogy of dust particles is similar to that in Volz (1972). In our study, we adopted the IR-to- 

20 VIS AOT ratio estimated from using the refractive index in Volz (1972) also because our 

21 focused domain is the eastern tropical North Atlantic Ocean. For single scattering albedo in the 

22 shortwave bands (ojo sw ), Forster et al. (2007) documented that the SSA for dust reported from 

23 various studies ranges from 0.90 to 0.99, with a median value of 0.96. In Myhre et al. (2003), 
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1 the observed SSA near Cape Verde is about 0.96-0.97. In this study, we set the SSA = 0.96. 

2 Other radiative properties including longwave single scattering albedo (coo' w ) as well as 

3 asymmetry factor in shortwave (g sw ) and longwave (g lw ) bands are set to values within the range 

4 reported in Dufresne et al. (2002) and Haywood et al. (2005). The complete list of settings for 

5 the prescribed dust radiative properties is given in Table 3. 

6 

7 4. Numerical simulations for 60 AEWs 

8 Using a stand-alone NASA/GSFC radiative transfer model (Chou and Suarez, 1999; 

9 2001), we found that the maximum radiative forcing exerted by the prescribed dust layer occurs 

10 at local noon. The instantaneous shortwave aerosol radiative forcing for both WAl and WAh is 

11 about -140 W m" at surface and -82 W m’“ at top-of-atmosphere (TOA). The longwave 

12 radiative forcing at surface is 28 W m’ 2 for WAl and 26 W m 2 for WAh; at TOA, 8 W m 2 for 

13 WAl and 11 W m ’ for WAh. As shown in Figure 4a, when the dust layer is placed between 700 

14 hPa and 850 hPa, the atmospheric temperature profile can be modified. The decrease in the net 

15 solar radiative flux, due to scattering, causes a net cooling of the air column below the dust layer 

16 by about 0.3 K day" 1 . Within the dust layer, the reduction in insolation is overtaken by 

17 absorption of solar radiation within the layer, resulting in a 3 K day’ 1 heating. On the other hand, 

18 dust longwave radiative forcing traps the terrestrial radiation and causes a warming of the air 

19 column below the dust layer by about 1.4 K day’ 1 . Within the dust layer, the dust aerosols cause 

20 a cooling effect due to longwave emission. As a result, at local noon, the prescribed dust layer 

21 causes a net 2 K day’ 1 warming in the dust layer, and about a 1 K day’ 1 warming for the air 

22 column below the layer. In Figure 4b, the dust-induced heating/cooling in the WAh 

23 configuration is slightly weaker and extends to higher altitudes than that the WAl configuration, 
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1 but the main characteristics are the same. Note that daily averaged dust-induced heating rate can 

2 be much smaller than the instantaneous dust-induced heating rate at local noon. The heating rate 

3 computed from our model is about the same magnitude to the heating rate computed in Carlson 

4 and Benjamin (1980), in which the authors reported that the daily mean shortwave clear-sky dust 

5 radiative heating rate within the dust layer for AOT = 1.0 is about 2 K day" 1 . Instantaneous 

6 clear-sky dust heating rate can be as large as 5 K day 1 as reported in Mallet et al. (2009) and 

7 Lemaitre et al. (2010), but long-term dust-induced heating rate climatology can be much smaller 

8 (Evan and Mukopadhyay, 2010). In addition, as aforementioned, uncertainties associated with 

9 dust radiative properties can result in different heating rate profiles. For example, Wong et al. 

10 (2009) reported that their heating rate is smaller than Carlson and Benjamin (1980), possibly due 

11 to different mineralogy and size distribution assumed in the model. 

12 As a result of the aerosol radiative forcing, the atmospheric temperature profile in the 

13 model can be significantly altered, as shown in Figure 5. Comparing the WAl and NA 

14 experiments (Figure 5a), the dust-induced heating causes the lower troposphere to warm by 

15 about 0.3-0.5 K below the dust layer and about the 0.2-0.4 K cooling within the dust layer, with 

16 an appreciable diurnal cycle. WAh experiments show a similar but weaker effect, about 0.2-0. 4 

17 warming below the dust layer and 0. 1-0.3 cooling within the dust layer (Figure 5b). The diurnal 

18 cycle of the temperature difference is due to the fact that aerosol shortwave radiative forcing 

19 follows the diurnal cycle of the solar irradiance. The fact that the magnitude of the temperature 

20 difference in Figure 5 is much smaller than the magnitude of the aerosol radiative heating rate in 

21 Figure 4 suggests a much smaller dust-induced heating/cooling over a day, compared with 

22 instantaneous heating at local noon (Figure 4), and a possible rapid adjustment of the 

23 temperature field by convection and circulation. Moreover, as shown in Figure 6, the alteration 
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1 of the temperature profile decreases convective stability after 24 hours of simulation for model 

2 spin-up. The decrease of convective stability is attributed to the warming of the lower 

3 troposphere. The magnitude of the tropospheric stability difference ranges from -0.1 K to -0.6 K 

4 with a diurnal cycle corresponding to that of the temperature difference. The change of 

5 temperature in the lower troposphere provides a favorable condition for shallow convection, 

6 which facilitates vertical mixing in the lower troposphere. As a result, the convective available 

7 potential energy (CAPE) over the domain increases by about 20%, from an average of 824 J kg’ 1 

8 for all the NA experiments to 1002 J kg 1 for the WAl and 1009 J kg 1 for WAh experiments. 

9 Hence, the isentropic surface is elevated in the WA experiments, and that causes the lower- 

10 tropospheric isentropic potential vorticity to increase for most cases (Figure 7). It is worth 

1 1 noting that the increase of CAPE in this study is the result of the radiative effect of the prescribed 

12 dust layer. While the thermal structure of the SAL, which is already embedded in the initial and 

13 boundary conditions in the simulations, can decrease CAPE and prohibit convections (Dunion, 

14 201 1), our results suggest that the dust-induced change of atmospheric temperature profile can 

15 counteract the effect of the SAL on modifying CAPE. 

16 To further understand the physical mechanism, we compute the Lorenz energy cycle 

17 (Lorenz 1955; 1967), which was originally developed to analyze the energetics of the global 

18 atmospheric circulation but has also been applied to analyze AEW disturbances (e.g., 

19 Karyampudi and Carlson, 1988; Hsieh and Cook 2007). Following the formulations in Hsieh 

20 and Cook (2007), we compute for the total atmospheric column (from surface to 100 hPa) the 

21 zonal available potential energy (P z ), zonal kinetic energy (K z ), eddy available potential energy 

22 (Pe), and eddy kinetic energy (Kr). We also calculate the energy conversion rate between P z and 

23 Pe (Ca), energy conversion rate between K z and K E (Ck), energy conversion rate between P E and 
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1 Ke (Cpk), and energy conversion rate between Pz and Kz (C z ). Finally, we compute the effective 

2 diabatic heating that produces Pz due to condensational heating (Hz ), the effective diabatic 

3 heating that produces Pz due to radiative heating (Hz R ), the effective diabatic heating that 

D 

4 produces Pz due to PBL scheme heating (Hz ), the effective diabatic heating that produces Pe 

5 due to condensational heating (He ), the effective diabatic heating that produces Pe due to 

6 radiative heating (H E R ), and the effective diabatic heating that produces Pe due to PBL scheme 

D 

7 heating (He ). For this calculation, we took the domain-average and discard the data from the 

8 first day of the simulations when the model is in the spin-up process and average the results 

9 between the 2 nd and last day of the simulation for all NA and WA simulations. 

10 In the 60 sets of simulations, we observe that the occurrence of stronger storms increases 

1 1 (Figure 8). The frequency of occurring the weakest AEWs, Ke less than 100 kJ, decreases 35 % 

12 in the WAl and 25% in WAh experiments, while the occurrence of the strongest AEWs remains 

13 the same. In Figure 9, it appears that the dust direct radiative forcing in WAl enhances 57 out of 

14 60 AEWs, ranging from a few percent to about 30 percent. In WAh experiments, 58 out of 60 

15 AEWs are intensified. However, the intensification in the WAh experiments is weaker than the 

16 WAl experiments. Furthermore, the AEWs intensify more rapidly over time (Figure not shown) 

17 as the dust radiative effect continuously builds up. 

18 The Lorenz energy cycle also reveals the physical processes involved in the 

19 intensification of AEWs. As shown in Table 4, except for Ca which slightly decreases, the other 

20 3 energy conversion rates increase. In particular, the 5.6% (5.1%) increase of eddy kinetic 

21 energy in WAl (WAh) experiments is directly attributed to, first and foremost, the enhanced 

22 conversion between the eddy available potential energy and the eddy kinetic energy, which 

23 increases 0.09 (0.063) W m' , about a 27% (19%) increase. The energy conversion between the 
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1 zonal kinetic energy and the eddy kinetic energy in WAl (WAh) experiments increases about 

2 0.01 (0.016) W m ". Furthermore, the diabatic heating terms suggest that the condensational 

3 heating is the primary energy source. The radiative heating plays a minor role and has a negative 

4 effect on the energy terms. This can be explained by the fact that as the storm gets intensified, 

5 the cloudy area expands and thus the cloud-radiation feedback causes a cooling effect on the 

6 atmosphere. As shown in Table 4, we think that the dust direct radiative forcing, albeit small, 

7 serves as a catalyst that destabilizes the atmosphere and facilitates convections. The increased 

8 cloud cover associated with the enhanced convections causes a negative feedback through 

9 radiation, but it is very small compared with the additional condensational heating released from 

10 the convections. The additional condensational heating increases available potential energy for 

1 1 both zonal mean flow and the eddy. In turn, the available potential energy is effectively 

12 converted to elevate the kinetic energy of the zonal mean flow and the eddy. The kinetic energy 

13 of the zonal mean flow can then be converted to increase the eddy kinetic energy. Among the 4 

14 energy conversion terms, we find that the local convective overturning plays the primary role to 

15 increase the eddy kinetic energy by converting the eddy available potential energy. 

16 It is interesting that the energy intensification due to the same aerosol radiative forcing 

17 varies and depends upon the particular character of each AEW. To determine which factors play 

18 a role, we examined the correlation of the intensification (AKe = Ke WA_ Kr NA ) with respect to 

19 environmental conditions typically considered as important factors for cyclogensis, such as SST, 

20 600 hPa relative humidity, CAPE, PBL height (Hpbl), vertical wind shear between 200 and 850 

21 hPa ( V shr ), total column precipitable water (Q co iumn), vertical temperature gradient, the initial 

22 strength of the eddy, etc. Among these variables, we found that the increase of eddy kinetic 

23 energy is significantly correlated with precipitable water (R = 0.50), vertical wind shear between 
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1 200 mb and 850 mb (R = 0.51), and PBL height (R = 0.49), for the WAl configuration. The 

2 multivariate regression model has the form: 

3 AKe= 0.05Qcoiumn + 0.97 V sh + 0.03H PBL - 1 15.8 , 

4 with the R-squared value equals 0.46. Other variables aforementioned are also moderately 

5 correlated with AKe with a correlation coefficient about or larger than 0.4, but they do not pass 

6 Student’s t-test with a significance level = 0.05 on top of the 3 variables. For WA H 

7 configuration, we found that the increase of eddy kinetic energy is significantly correlated with 

8 precipitable water (R = 0.49) and PBL height (R = 0.46), and the multivariate regression model 

9 has the form: 

10 AKe= 0.05Qcolumn + 0.03HpBL — 107.6 , 

1 1 with the R-squared value equals 0.36. It should be noted that the development of AEWs depends 

12 on the environmental conditions as well as the interaction between the AEW and the 

13 environment, which can vary case by case. The dust-induced intensification in this study is 

14 found to be associated with enhanced local convective overturning. Thus, when the dust 

15 radiative forcing fortifies the original environmental conditions to facilitate local convective 

16 overturning, the dust-induced intensification is found the most evident. From the multivariate 

17 regression model, it is suggested that when the environment is moist and has an unstable PBL, 

18 which contribute to convective instability, the dust radiative forcing is likely to intensify the 

19 AEW more effectively. For the WAl experiments, the larger baroclinicity, which is associated 

20 with the north-to-south temperature gradient, can also contribute to larger intensification of 

21 AEWs since the dust-induced heating further enhances the horizontal temperature gradient. In 

22 WAh experiments, vertical wind shear was also well- correlated with the increase of eddy kinetic 
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1 energy (R=0.41), but it did not pass the significant test to be included in the regression model on 

2 top of precipitable water and PBL height. 

3 

4 5. Conclusions 

5 In this study, we performed a large number of numerical simulations to explore the 

6 possible role of radiative forcing by mineral dust in intensifying AEWs. While acknowledging 

7 the large uncertainty associated with the radiative properties and distribution of dust aerosols, we 

8 applied prescribed dust layers in the model with radiative properties and distributions within the 

9 range reported in previous studies. The simplified dust layer helps us to isolate the direct 

10 radiative effect of dust on the AEW development from non-linear responses of the aerosol- 

1 1 meteorology interaction in a model that fully couples the aerosols, clouds, radiative transfer, and 

12 circulation. 

13 Despite the uncertainties associated with mineral dust, the large number of simulations 

14 using prescribed simple dust layers establishes the statistical validity of the radiative effect of the 

15 dust layer on the intensification of AEWs. We found that dust radiative forcing modifies the 

16 vertical temperature profile by warming the lower tropopsphere and thereby decreases 

17 atmospheric stability. As stability decreases, most AEWs intensify. The model results are 

18 consistent with the observational studies of Alpert et al. (1998) and Jones et al. (2004). 

19 Furthermore, the Lorenz energy cycle analysis suggests that in the presence of a dust layer, the 

20 condensational heating increases dramatically and is effectively converted to eddy kinetic energy 

21 through convective overturning. The magnitude of intensification varies from one AEW to 

22 another and is the strongest when the environment is characterized as moist, baroclinic, and has a 

23 higher PBL. 
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5 (middle column); and (3) Temperature (color shaded, unit: K) and winds (vector, unit: m s' 1 ) at 
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9 Figure 2. Time series of ECMWF relative vorticity at 700 hPa (purple striped area encloses 

10 relative vorticity larger than 1.0 x 10 5 s' 1 ) and MODIS AOT (yellow shaded area), zonally 

11 averaged between 15°W and 50°W, in (a) 2003; (b) 2004; (c) 2005; (d) 2006; and (e) 2007. Red 

12 stars denote the dates and latitudes of TC-genesis identified by National Hurricane Center. 
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14 Figure 3. Composite of MODIS AOT (color shaded), and ECMWF relative vorticity (contour; 

15 interval = 0.5 x 10" 5 s' 1 ) and winds (vector; unit: m s' 1 ) at 700 hPa for the 60 AEWs on (a) the 
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19 Figure 4. Aerosol-induced clear-sky instantaneous heating rate (unit: K/day) at 12:00 (local 

20 time). Dotted area denotes the location of the prescribed dust layer. 
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1 Figure 5. Domain- averaged time-height plot of temperature (K) difference between (a) WAl and 

2 NA experiments; and (b) WAh and NA experiments, averaged over 60 pairs of simulations. 

3 Dotted area denotes the location of the prescribed dust layer. 
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5 Figure 6. Time series of domain-averaged tropospheric stability difference (defined as potential 

6 temperature difference between 200 hPa and 1000 hPa) between 60 pairs of (a) WAl and NA 

7 experiments; and (b) WAh and NA experiments (unit: K). 
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9 Figure 7. Time series of domain- averaged 315 K potential vorticity difference between 60 pairs 

10 of (a) WAl and NA experiments; and (b) WAh and NA experiments (unit: PVU). 
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12 Figure 8. Histogram of eddy kinetic energy (unit: kJ/m') of 60 AEWs in (a) NA (solid fill) and 

13 WAl (crosshatch) experiments, and (b) NA (solid fill) and WAh (crosshatch) experiments. 
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15 Figure 9. Histogram of eddy kinetic energy increase (unit: percent) of the 60 AEWs due to 

16 aerosol radiative forcing at 700-850hPa (solid fill) and 600-850hPa (crosshatch). 
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AEW disturbance 

Dust outbreak 
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2004 
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2005 

32 
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2006 
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3 
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32 
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2 Table 1. Number of AEW disturbances, dust outbreaks, TC-genesis events in 2003, 2004, 2005, 

3 2006, and 2007. 
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2 Table 2. List of the starting date of AEW cases. 
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sw 
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„ sw 

COo 

gSW 

Iw 
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Iw 
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1.0 

0.96 

0.7 

0.5 

0.2 

0.4 


Table 3. Aerosol radiative properties including optical depth in shortwave (r sw ), single scattering 
albedo in shotwave (coo sw ), asymmetry factor in shortwave (g sw ), optical depth in longwave (x lw ), 
single scattering albedo in longwave (coo lw ), and asymmetry factor in longwave (g lw ). 
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Simulation 

NA 

WA l 

WA h 

Energy terms 
(kJ m 2 ) 

Pz 

30.5 

27.3 

29.1 

Kz 

183.3 

182.2 

185.6 

Pe 

35.8 

35.4 

35.2 

K e 

126.6 

133.9 

133.0 

Energy conversion rates 
(W m 2 ) 

C z 

-0.081 

-0.024 

-0.031 

C A 

-0.007 

-0.008 

-0.008 

C k 

0.068 

0.078 

0.084 

CpK 

0.332 

0.422 

0.395 

Diabatic heating terms 
(W m' 2 ) 

ESI 

-0.012 

0.046 

0.015 

ESI 

-0.044 

-0.082 

-0.072 

ESI 

-0.008 

-0.009 

-0.007 

EM 

0.197 

0.256 

0.241 

ESI 

0.090 

0.033 

0.043 

ESI 

-0.005 

-0.005 

-0.005 


1 

2 Table 4. Energy terms, energy conversion rates, and diabatic heating terms averaged over 60 

3 AEWs. 
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2 Figure 1. Five-month (June-October) average (1) MODIS AOT (yellow shaded), and positive 

3 relative vorticity (contoured, unit: 10 5 s' 1 ) and winds (vector, unit: m s' 1 ) at 700 hPa (left 

4 column); (2) Temperature (color shaded, unit: K) and winds (vector, unit: m s' 1 ) at 700 hPa 

5 (middle column); and (3) Temperature (color shaded, unit: K) and winds (vector, unit: m s' 1 ) at 

6 1000 hPa (right column) for years 2003 (top)-2007 (bottom), respectively. Relative vorticity, 

7 winds, and temperature are obtained from ECMWF operational analysis. 
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Aerosol Optical Thickness Relative Vorticity TC-genesis 



I I I 1—1 

0.4 0.6 0.8 1 


2 Figure 2. Time series of ECMWF relative vorticity at 700 hPa (purple striped area encloses 

3 relative vorticity larger than 1.0 x 10 2 3 4 5 s' 1 ) and MODIS AOT (yellow shaded area), zonally 

4 averaged between 15°W and 50°W, in (a) 2003; (b) 2004; (c) 2005; (d) 2006; and (e) 2007. Red 

5 stars denote the dates and latitudes of TC-genesis identified by National Flurricane Center. 
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(a) Day 1 
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Figure 3. Composite of MODIS AOT (color shaded), and ECMWF relative vorticity (contour; 
interval = 0.5 x 10 5 s' 1 ) and winds (vector; unit: m s' 1 ) at 700 hPa for the 60 AEWs on (a) the 
day of; (b) 1 day after; (c) 2 days after; and (d) 3 days after, the AEW s leave the coast of African 
continent. The thick purple line denotes zero relative vorticity. 
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Clear-sky instantaneous heating rate 


(a) AOT sw =1 .0 700-850hPa 



1 

2 Figure 4. Aerosol-induced clear-sky instantaneous heating rate (unit: K/day) at 12:00 (local 

3 time). Dotted area denotes the location of the prescribed dust layer. 
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1 

2 Figure 5. Domain- averaged time-height plot of temperature (K) difference between (a) WAl and 

3 NA experiments; and (b) WAh and NA experiments, averaged over 60 pairs of simulations. 

4 Dotted area denotes the location of the prescribed dust layer. 
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Tropospheric Stability Difference 
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1 

2 Figure 6. Time series of domain-averaged tropospheric stability difference (defined as potential 

3 temperature difference between 200 hPa and 1000 hPa) between 60 pairs of (a) WAl and NA 

4 experiments; and (b) WAh and NA experiments (unit: K). 
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PV 315 k Difference 
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1 

2 Figure 7. Time series of domain- averaged 315 K potential vorticity difference between 60 pairs 

3 of (a) WAl and NA experiments; and (b) WAh and NA experiments (unit: PVU). 
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Eddy Kinetic Energy 

(a) 700-850hPa 
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2 

2 Figure 8. Histogram of eddy kinetic energy (unit: kJ/nr) of 60 AEWs in (a) NA (solid fill) and 

3 WAl (crosshatch) experiments, and (b) NA (solid fill) and WAh (crosshatch) experiments. 
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Eddy Kinetic Energy Increase 



1 

2 Figure 9. Histogram of eddy kinetic energy increase (unit: percent) of the 60 AEWs due to 

3 aerosol radiative forcing at 700-850hPa (solid fill) and 600-850hPa (crosshatch). 
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